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We study the ^So proton pairing gap in neutron star matter putting emphasis on influence of the 
Dirac effective mass and the proton fraction on the gap within the relativistic Hartree-Bogoliubov 
model. The gap equation is solved using the Bonn-_B potential as a particle-particle channel inter- 
action. It is found that the maximal pairing gap Amax is 1-2 MeV, which has a strong correlation 
with the Dirac effective mass. Hence we suggest that it serves as a guide to narrow down parameter 
sets of the relativistic effective field theory. Furthermore, the more slowly protons increase with 
density in the core region of neutron stars, the wider the superfiuid range and the slightly lower the 
peak of the gap become. 
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I. INTRODUCTION 

Superfluidity in neutron star matter is one of the hot 
issues in nuclear astrophysics since superfluidity plays a 
key role in affecting the cooling of neutron stars 0] ■ Un- 
like the neutron pairing, the proton pairing oc- 
curs in dense matter with supranuclear density, where its 
properties are not much known. Such a region is highly 
relevant to the URCA processes that control cooling of 
neutron stars. 

In neutron stars, several types of baryon pairing are be- 
lieved to appear. In the inner crust region, neutrons will 
form the ^5*0 pairs j2|,|^|^]. At the corresponding baryon 
density 10~^po ^ Pb ^ 0-7po, where po — 0.15 fm~^ 
is the saturation density of symmetric nuclear matter, 
the ^5*0 partial wave of the nucleon-nucleon (NN) inter- 
action is attractive. This attraction helps neutrons pair 
up through the well-known BCS mechanism. In the core 
region ps > O.Tpo, the '^P2 neutron pairing may also ap- 
pear since the '^P2 partial wave of the NN interaction 
becomes attractive enough 0, 13- In contrast, the ^Sq 
partial wave would become repulsive there so that the 
neutrons would cease to pair in the ^Sq state. Instead, 
the ^Sq proton pairs are predicted to appear owing to 
its fraction smaller than neutrons nil. At much more 
higher baryon density pB ^ Spo, various hyperons may 
emerge; some kinds of them possibly form pairs in the 
same way as nucleons do @, 0i 0- We, however, stay 
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on a picture of neutron stars without hyperons in the 
present study. 

The size of the nuclear pairing gap in dense matter is 
controversial since many uncertainties remain regarding 
NN interactions in dense medium, methods of approxi- 
mation, sparsity of experimental data under extreme con- 
ditions, and so on. From the perspective of the present 
status of immaturity, a great deal of study with diversity 
is absolutely necessary. In the meantime, many studies of 
neutron stars have been performed using various frame- 
works so far. Recently, relativistic models are attracting 
attention in researches on neutron stars since they are 
suitable to describe the stars in compliance with the spe- 
cial relativity 0. Most often used among them is the 
relativistic mean field (RMF) model, particularly owing 
to its economical way of description. Hence we choose the 
extended model of it, the relativistic Hartree-Bogoliubov 
(RHB) model as the framework of the present study. 

The primary aim of this paper is to elucidate effects 
of bulk properties on the ^5*0 proton pairing correlation 
in neutron star matter (consisting of n, p, e~ , and p,~) 
using the RHB model with capability of handling the 
pairing correlation. Here we show importance of environ- 
mental properties of dense matter that surrounds Cooper 
pairs of protons (the pairs themselves are of course a part 
of the environment). In particular, the Dirac effective 
mass of nucleons and the proton fraction are considered 
as accompanying quantities of great importance since the 
bulk properties are influential on superfluidity in neutron 
stars. Therefore we aim to address the bulk properties 
of neutron star matter and its superfluidity on the same 
footing within the RHB model. 

The present study covers the two aforementioned quan- 
tities: The Dirac effective mass has an effect on the pair- 
ing correlation via the density of states. The symmetry 
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energy coefficient controls the proton fraction in neutron 
star matter, especially its density dependence. Needless 
to say, they connect with the equation of state (EOS) of 
dense matter; it is strongly related to properties of neu- 
tron stars, such as their internal structure, mass, radius, 
and so on. These macroscopic properties of neutron stars 
set stringent microscopic requirements for constituents; 
hence the pairing correlation is no exception. In this 
study, we use two distinct RMF models with two dis- 
tinct parameter sets for each to compose neutron star 
matter; thereby we have four distinct EOSs here. Two 
of them are for a comparison of the effect of the Dirac 
effective mass, the other two are for a comparison of the 
effect of the proton fraction. 

This paper is organized as follows. In Sec. model 
Lagrangians and the gap equation for the ^5*0 proton 
superfluidity are illustrated. In Sec. IIIII results of the 
pairing properties in neutron star matter are presented. 
Section ITVI contains our summary and conclusions. 



II. MODELS 

As is well known nowadays, the RMF model La- 
grangian with nonlinear cubic, quartic terms of ct bosons, 
and a quartic one of u mesons has achieved remarkable 
successes in studies on nuclear/hadronic physics. This 
type of Lagrangian was initiated by .Bodmer 10] and has 
matured steadily 0, . A parameter set is determined 
so as to reproduce the saturation properties of symmet- 
ric nuclear matter and ground-state properties of typical 
finite nuclei. A prominent feature of the model is that 
scalar and vector self-energies (or Hartree fields) follow 
the density-dependence of those obtained by the Dirac- 
Brueckner-Hartree-Fock (DBHF) approach even in mod- 
erately supranuclear density. Hence, the model is con- 
sidered to have the capability to deal with neutron star 
matter at such density. This model Lagrangian is of the 
form 



+ i(a,a)(aM-^->^-f-^-f-^ 



(1) 



where fi^i, = d^j^uju ~ d^uj^ and B^^ = d^b„ - 9^6^. 
The symbols -0, cr, lo^, h^, M, mo-, m^^, and rup signify 
the fields of nucleons, a bosons, lo mesons, p mesons, 
the masses of nucleons, a bosons, oj mesons, and p 
mesons, respectively. We call Eq. "standard RMF 
Lagrangian." 

As an extension of the well-established standard RMF 
model, the effective field theory (EFT) provides a 
modern aspect of the RMF model; an energy density 



functional obtained from the RMF Lagrangian approxi- 
mates the exact energy density functional of the ground 
state of a hadronic system in the sense of the density func- 
tional theory (DFT). Adding apt interaction terms, ab- 
sent in the standard one, to the RMF functional advances 
it to the exact functional little by little. In this EFT- 
inspired approach, mesons and baryons are not taken 
as elementary degrees of freedom, so that a Lagrangian 
composed of corresponding fields is allowed to be non- 
renormalizable. This leads to, in principle, unrestricted 
inclusion of any terms of meson sclf-intcractions consis- 
tent with symmetries of the underlying theory, QCD. 
These terms are completely absent in the standard RMF 
Lagrangians that respect renormalizability. Recently, 
iHorowitz and Piekarewic^ proposed an RMF Lagrangian 
along with concepts and methods of the EFT and the 
DFT. They extended the standard RMF Lagrangian with 
a nonlinear uo-p coupling particularly from the viewpoint 
of the symmetry energy. The new coupling paves the 
way to change the density dependence of the symme- 
try energy, and hence the proton fraction in neutron 
star matter. For the extended Lagrangian, the present 
study employs an additional interaction Lagrangian of 
the form Q 



/Ieft = C + Ag^bp ■ b^" A^gluj^oj" , 



(2) 



where A^ signifies a nonlinear coupling constant. We 
name Eq. |5J "EFT-inspired Lagrangian." The impor- 
tant thing is that the Lagrangian ||2J) gives the symmetry 
energy coefficient 
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The symbol m* denotes the effective mass of p mesons 

mf^ml + 8glA^gl{oj„f, (4) 

where {O) represents an expectation value of a field O 
for the ground state of the system. In this study, we 
pick out the NL3-hp and Z271 parameter sets [l4 [l5|: 
each set gives the lower or upper limit of commonly ac- 
cepted range of the effective nucleon mass in symmetric 
nuclear matter at po, M* ~ (0.6-0.8)Af. Most RMF pa- 
rameter sets yield the mass within this range. Note th at 
we rename the NL3 set by iHorowitz and Piekarewic^ in 
Ref. here because of their slight modification of the 
well-known NL3 parameter s et jlq . We will also men- 
tion the TMl parameter set |l2j for the standard RMF 
Lagrangian in addition to the above two sets. Table HI 
shows the parameter sets. 

In neutron star matter composed by the models, we 
solve the gap equation for the ^5*0 proton pairing correla- 
tion at zero temperature using the Bonn-i? potential [T^ 
as a particle-particle (p-p) channel interaction w, 



A(p) = - 



I 

8^ 



A(fc) 



X v{p, k) k^dk, (5) 
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TABLE I: Parameter sets used in the present study. We fix the bare nucleon mass M = 939.0 MeV and the p meson mass 
m,p = 763.0 MeV, except for TMl where M = 938.0 MeV and m,p = 770.0 MeV are used. 



Parameter set 


m„ [MeV] 


m„ [MeV] 


9CT 




92 [fm-i] 


53 


C3 


9p 


NL3-hp (A^ = 0) 


508.194 


782.5 


10.217 


12.868 


-10.431 


-28.885 





4.461 


NL3-hp (A^ = 0.025) 


508.194 


782.5 


10.217 


12.868 


-10.431 


-28.885 





5.376 


Z271 (A„ = 0) 


505.0 


783.0 


7.031 


8.4065 


-5.4345 


-63.691 


49.94 


4.749 


Z271 (A„ = 0.040) 


505.0 


783.0 


7.031 


8.4065 


-5.4345 


-63.691 


49.94 


5.008 


TMl (A^ = 0) 


511.198 


783.0 


10.0289 


12.6139 


-7.2325 


0.6183 


71.3075 


4.6322 



where Ek = y/WTu*^ + gu,{uJa) + gMo'^)- The p-p 
channel interaction between protons v{p,k) is obtained 
through angular integration with respect to the angle be- 
tween linear momenta p and k. This serves as a process 
of projecting out the S-wave component of the interac- 
tion. Its original form is nothing but the antisymmetrized 
matrix element of the employed interaction V , which is 
defined by 

5(p, k) = (ps', ps'|F|ks, ks) — (ps'j ps'|V^|ks, ks), (6) 

where a tilde over its argument denotes time reversal. 
We use the lowest approximation as to the p-p chan- 
nel interaction to study exclusively the effect of the bulk 
properties on superfluidity. It is, however, known that 
the polarization effects bring about significant reduction 
of the gap within nonrelativistic frameworks 0, . 

The models employed here should be regarded as hy- 
brid models because we use different interactions in the 
particle-hole {p-h) and the particle-particle channel. This 
prevents ambiguity in choosing the pairing interaction. 
Thus we concentrate on the study of the effect of the 
Dirac effective mass and the proton fraction on the pro- 
ton superfluidity. 



III. RESULTS AND DISCUSSION 

To begin with, we calculate the ^5*0 pairing gaps and 
the Dirac effective masses of protons without the isovec- 
tor nonlinear coupling. The results are shown in Fig. ^ 
which apparently shows that the gaps strongly depend 
on the Dirac effective masses. The smaller the mass is, 
the smaller the pairing gap is. The NL3-hp set makes 
symmetric nuclear matter saturated with the Dirac ef- 
fective mass of 0.59M, while the Z271 set does with that 
of 0.80M. According to Ref. both sets give similar 
properties of typical finite nuclei and low-mass neutron 
stars on the one hand, different properties of typical neu- 
tron stars of mass 1.4A/0 follow from the sets on the 
other hand. For comparison, we also present the result 
of the TMl set in Fig. with the light gray curves. The 
effective mass obtained by the TMl set lies in between 
those by the other two sets, so does the pairing gap; it is 
consistent with the above statement as well. In view of 
the pairing correlation in neutron star matter, using the 
different p-h interactions and the same p-p interaction 
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FIG. 1: ^Sq proton pairing gaps (left scale, solid curves) and 
Dirac effective masses (right scale, dashed curves) as functions 
of baryon density in neutron star matter using the standard 
RHB model with the NL3-hp, Z271 [3, and TMl [T^ pa- 
rameter sets. 



gives the significant difference shown in Fig. ^ Hence, 
this suggests that we can use the strong correlation as a 
guide to narrow down the parameter sets favorable for a 
calculation such as the present study. 

Next, we show the proton fractions Yp = Pp/ Pb in 
Fig. 12 for the variations of Lagrangian. The solid black 
curve denotes the proton fraction obtained by the NL3-hp 
set, the solid gray curve by the Z271 set, both without 
the new nonlinear coupling. Their characteristic is the 
large proton fraction, which is the typical result obtained 
from the standard RMF Lagrangian. The result obtained 
by the TMl set is similar to the one by the NL3-hp set 
without the coupling (and hence is omitted hereafter). In 
contrast, using finite A„, namely, the EFT- inspired RMF 
Lagrangian, we have obtained smaller proton fractions 
at densities corresponding to the core region of neutron 
stars, Pb ^0.1 fm^'', than using A^^ = 0; these fractions 
are drawn with dashed curves in Fig. |21 As the original 
intent of this kind of Lagrangian 14], it yields gentle in- 
crease of the fraction at these relevant densities owing to 
the symmetry energy adjustable via the isovector nonlin- 
ear coupling indicated by Eqs. © and Q. Not shown 
are curves of the proton fraction that lie in between the 
two extreme conditions A^^ = and 0.025 (0.040) for the 
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FIG. 2: Proton fractions of the four distinct models of neutron 
star matter as functions of baryon density. The thin region 
hatched in gray represents a threshold of the proton fraction 
for the direct URCA process, about 13-15 %. The number in 
parenthesis in the legend indicates the value of A^^. 




FIG. 3: ^5*0 proton pairing gaps as functions of the proton 
Fermi momentum in neutron star matter using the RHB mod- 
els with the NL3-hp and Z271 parameter sets with and with- 
out Kuj . The legend is the same as in Fig. |21 



NL3-hp (Z271) set, when it is varied between them. To 
wind up the account of Fig. [3 we should note how each 
curve of the proton fraction depends on baryon density, 
with its influence on the pairing correlation in mind. 

To offer more direct information on the ^5*0 proton 
pairing gaps at the proton Fermi surface and their re- 
lation to the proton fractions, we present the gaps as 
functions of the proton Fermi momentum in Fig. |31 An 
apparent difference between the gaps obtained with the 
NL3-hp and Z271 sets stems, again, from the difference of 
the Dirac effective masses between the sets, as shown in 
Fig.n III addition, the two curves for each parameter set 
with and without A^ reflect the increase rate of the pro- 
ton fraction as functions of baryon density (see Fig. [SJ . 
Comparing the property of protons at the same Fermi 
momentum explains further details about that: With A^, 
where the rate is moderate at the relevant densities, the 
proton Cooper pairs are immersed in denser background 
than those without Ai^, where the rate is rapid. Mean- 
while, the higher the baryon density is, the smaller the 
Dirac effective mass of protons is. We thus obtain the 
smaller pairing gaps at the same Fermi surface for the 
parameter sets with A^^, namely, for the EFT-inspired 
Lagrangian. 

We compare the pairing gaps as functions of baryon 
density obtained with and without the isovector nonlin- 
ear coupling Aij in Fig. 0] Since the density dependence 
of the Dirac effective masses obtained from the respec- 
tive parameter sets is almost the same for any values of 
A(^, the obtained gaps reflect the respective proton frac- 
tions presented in Fig. 13 As a whole, the pairing gaps 
survive in higher density region for the EFT-inspired La- 
grangian with flnite A(^, and have lower peaks than for 
the standard Lagrangian with Aj^ = 0. Figures |21 and ^ 
give indications as follows. 
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FIG. 4: ^ Sii proton pairing gaps as functions of baryon den- 
sity in neutron star matter using the RHB models with the 
NL3-hp and Z271 parameter sets with and without . Solid 
curves and the legend are the same as in Figs. and |21 re- 
spectively. 

(1) Up to pb — 0.05 fm^"^, the proton fractions are 
almost the same except the one of the NL3-hp set with 
flnite A^, which is a little larger than the others (see 
Fig. [2Jl . The enhancement of the gap due to the large 
proton fraction cancels out the smallness of the effective 
mass given by the NL3-hp set. Hence, the pairing gap 
obtained by the NL3-hp set with A(^ — 0.025 is almost 
the same as those by the Z271 set. We note, however, 
that this region roughly corresponds to the inner crust 
of neutron stars, where the ^5*0 proton pairing is much 
weaker than the neutron counterpart. 

(2) At Pb ~ 0.1 fm~^, as the proton fractions of the 
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FIG. 5: Critical temperatures of ^So proton superfluids as 
functions of baryon density in neutron star matter using the 
RHB models with the NL3-hp and Z271 parameter sets with 
and without A^^ . The legend is the same as in Fig. |5| 



Z271 set start to deviate from each other, so do the pair- 
ing gaps obtained by the set. As to the NL3-hp set, the 
proton fractions for A^^ = and 0.025 nearly coincide 
there, so do the pairing gaps. 

(3) For 0.1 fm-3 < ps < 0.2 fm^^, the high pro- 
ton fractions make the corresponding pairing gaps large. 
This implies that richness of protons favors the pairing 
correlation by taking advantage of an attractive part of 
the proton-proton interaction in this region. 

(4) At pb ~ 0.2 fm^^, an opposite situation to the pre- 
vious case (3) is realized; namely, the high proton frac- 
tions make the corresponding pairing gaps small. This 
implies that richness of protons disfavors the pairing cor- 
relation by virtue of a repulsive part of the proton-proton 
interaction in this region. As a result of the difference be- 
tween the effective masses, the pairing gap given by the 
Z271 set is still very large. To sum up the discussion 
of Fig. ^ the more gently protons multiply in the core 
region of neutron stars, up to the higher density super- 
fluidity survives and the somewhat smaller the peak of 
the gap one obtains. 

Let us here present three figures relevant to the physics 
of neutron stars. One is Fig.jS] which represents the den- 
sity dependence of the critical temperatures of the super- 
fluids. The critical temperatures Tc are obtained by the 
universal relation for a weak-coupling BCS superconduc- 
tor at T = K, fceTc = 0.57A(4''^ T = 0), where fee 
is the Boltzmann constant. Since the temperature inside 
evolved neutron stars is about 10^ K, Fig. |31 signifies that 
the superfluidity is likely to exist in the all cases consid- 
ered. The other two are Figs. El and [7| which depict the 
neutrino emissivities Q and their reduction factor R by 
the ^5*0 proton superfluidity as functions of baryon den- 
sity, drawn at the internal temperatures T = 1.0 x 10^ K 
and 3.0 x 10* K. We now take, for example, the nucleon 
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FIG. 6; Neutrino emissivities of the nucleon direct URCA 
process (left scale, solid curves) and the reduction factor for 
^So proton pairing (right scale, dashed curve) as functions 
of baryon density in neutron star matter for the NL3-hp pa- 
rameter set without A^^. Solid gray curves represent the neu- 
trino emissivities in normal fluid. The upper panel is drawn 
for the internal temperature T = 1.0 x 10* K, the lower for 
T = 3.0 X 10** K. 
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FIG. 7: Same as Fig.O but for the Z271 parameter set with- 
out K^. Also note that the right scales are different from 
those of Fig. |S1 



direct URCA process as a cooling agent of neutron stars. 
The neutrino emissivity of the nucleon direct URCA pro- 
cess in normal fluid Qo has the following form in units of 
erg s"^ cm~3 |23: 



)o = 7.55 X lO^'VeTf 



MfMf 



M2 



— 9{pYe+P¥p-P¥n), {^) 
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where is the chemical potential of electrons (and 
muons due to the condition of beta equilibrium), Tg is 
the temperature in units of 10^ K, and 9{pFe +PFp—PFn) 
is the triangle condition for particle momenta. After the 
appearance of muons, they also contribute to the neu- 
trino emission so that the cmissivity Qo just doubles, 
which manifests itself as the small kinks after the igni- 
tion of the process involving electrons (abrupt increase of 
the emissivities) in Figs. El and [7| It should be noted that 
we approximate the effective masses M* = M* = M* in 
Eq. ||7J|: To distinguish between them in the RMF model, 
we must introduce isovector-scalar mesons and shortly 
discuss this issue later. The neutrino emissivity in super- 
fluid is thus written as Q = QqR, the reduction factor 
used here being 



R 



exp 
exp 



-A(4P\T)/fcBT 



-A(4P';T = 0)/fcBT 



exp -1.76r,(4''^)/T , 



(8) 



which is the form often used in the past studies. It should 
be mentioned that the reduction factor of the form JSJ 
leads to an overestimate of the suppressive effect 21, 22]. 
Although we have employed the pairing gaps at T = 
K for a qualitative estimate in the present study, one 
should use the pairing gaps at finite temperatures for a 
quantitative study of thermal evolution of neutron stars. 

As shown in Figs. ^ and 13 both the NL3-hp and the 
Z271 set with A„ — yield the relatively large proton 
fraction. Meanwhile, the pairing gap for the NL3-hp set 
closes around the density at which the direct URCA pro- 
cess takes effect. Therefore, the reduction effect is not 
much large, which is clear for T = 3.0 x 10^ K as shown 
in Fig. IB] By comparison, the superfluid range for the 
Z271 set is wide enough to cover the density region where 
the direct URCA process is turned on. Hence the process 
is strongly suppressed by the interior proton superfluid 
as shown in Fig. [7| For both sets with finite A(^ , on the 
other hand, the density regions of the superfluidity and 
the direct URCA process do not overlap each other; the 
latter starts at ps = 0.442 fm'^ with Yp = 0.136 for the 
NL3-hp set and at ps = 1.028 fm'^ with Yp = 0.141 for 
the Z271 set. Thus the superfluidity is of no suppressive 
effect on the direct URCA process for both with finite 
A„. 

Now we develop a brief discussion on the cooling of 
neutron stars relative to the above results. From the as- 
pect of surface temperature of neutron stars, we can clas- 
sify observed neutron stars into hotter and colder ones. 
Broadly speaking, this means that there exist two major 
cooling scenarios; the so-called "standard" and "nonstan- 
dard" cooling j23| . The standard cooling is dominated by 
the modified URCA process, which cools neutron stars 
slowly and thus results in hotter stars. In contrast, the 
nonstandard cooling includes the nucleon/hyperon di- 
rect URCA processes, via which neutron stars cool faster 



and lead to colder stars, in conjunction with the sup- 
pression due to superfluidity. Concerning this classifica- 
tion, the pulsar PSR J0205-I-6449 was recently discov ered 
in the supe rnova remnant 3C 58 by 'M urrav et all [23| . 
ISlane et al\ subsequently deduced a strong upper limit 
on the surface temperature 25]. They showed that the 
limit falls below the prediction of the standard cooling 
model, which suggests that PSR J0205-)-6449 seems to 
follow the nonstandard cooling scenario. In the first 
place, the proton fraction higher than about 13-15 % 
(shown by the hatched region in Fig. O activates the 
direct URCA process involving nucleons in neutron star 
matter with the standard RMF Lagrangian. However, 
our results from the Z271 set show that the proton pair- 
ing gap is large enough to suppress the process, so that 
the resulting cooling of neutron stars is not much effec- 
tive to meet observational data of colder neutron stars. 
This calls for the cooling that results from hyperons or 
meson condensates; the latter consequently means coex- 
istence of superfiuidity and meson condensations. In the 
second place, using the EFT-inspired RMF Lagrangian 
with the Z271 set, neutron star matter is composed with 
the proton fraction lower than the direct URCA thresh- 
old up to ~ 7po- This prohibits neutron stars from cool- 
ing by the direct URCA process involving nucleons; the 
modified URCA process dominates their cooling instead, 
given the present constitu ent particles. In this regard, 
iHorowitz and Piekarewic j studied the direct URCA pro- 
cess in detail using the EFT-inspired RMF Lagrangian 
with intermediate values of A^^ |26j . Bear in mind that 
we have put aside the Cooper pair breaking and forma- 
tion processes |2^ that work as a cooling accelera- 
tor, while we have taken superfluidity into account as a 
cooling retardant alone: Since the neutrino emissivity of 
these processes sensitively depends on the density and 
the temperature dependence of the pairing gap, our ap- 
proximate treatment at the moment does not allow us 
to estimate this emissivity and to determine the scenario 
that the pulsar follows. 

Finally, we should note preceding works on the ^Sq 
proton superfluidity in neutron star matter. Accord- 
ing to them, nonrelativistic approaches and the Dirac- 
Brueckner-Hartree-Fock approach give similar results 
with the maximal pairing gap A^ax S, 1-0 MeV and the 
closing density pc ~ 0.4 fm'^ plIsoTlp]. The discrep- 
ancy between them and our results having been exposed, 
with the reservation of the lowest approximation, we now 
need information to narrow down further RMF param- 
eter sets for the purpose like the present study. In this 
context, extracting the EOS from direct observation of 
the stars and experiments of neutron-rich nuclei is indis- 
pensable. 



IV. SUMMARY AND CONCLUSIONS 

We have investigated the ^Sq proton pairing in neutron 
star matter using the relativistic Hartree-Bogoliubov 
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model. Since proton Cooper pairs reside in extremely 
dense surroundings inside neutron stars, special care 
should be taken of the effects of the bulk properties, as 
well as pairing interactions and levels of the approxima- 
tion. Therefore, we have studied the effects on the pairing 
correlation using the lowest approximation with respect 
to the pairing interaction. We summarize and conclude 
the following: First, using the standard RMF parame- 
ter, we have obtained that the maximal pairing gap is 
about 1-2 MeV, dependent on the parameter sets used, 
which is comparable or larger than the values obtained 
in the preceding studies. This clearly shows the strong 
correlation between the effective mass of nucleons and 
the pairing gap, which was also concluded within nonrel- 
ativistic models. Although the relativistic model tends 
to have a smaller effective mass of nucleons (the Dirac 
effective mass) than those in nonrelativistic models, the 
value of the maximal gap is not much small against gen- 
eral expectations. Second, we have studied the problem 
also using the extended L agrangian proposed in line with 
a concept of the EFT bv lHorowitz and Piekarewicd [l^ . 
The proton fraction obtained can be controlled by the ad- 
ditionally introduced parameter . We have found that 
the more slowly protons increase with density in the core 
region of neutron stars, the wider the superfluid range 
and the slightly lower the peak of the gap become. 

In this study, we have ruled out direct effects of the 
bulk properties on the p-p channel interaction: For in- 



stance, self-consistent (or in-medium) Dirac spinors in 
the p-p channel are not used here though they proba- 
bly affect properties of the pairing correlation at high 
density Q. Furthermore, we have included no isovcctor- 
scalar meson, which may play important roles in high- 
density and isovector physics through larger effective 
mass of protons than that of neutrons [s^l- In general, 
the isovector-scalar meson increases protons in the high- 
density region, which is an effect opposite to that of the 
EFT-inspired Lagrangian employed here. Last but not 
least, having clarified the effects of the bulk properties 
on superfluidity by solving the gap equation at zero tem- 
perature, we are prepared for solving it at finite tempera- 
tures. It will be mandatory to consider these issues. With 
attention to the bulk properties of neutron star matter, a 
study of baryon superfluidity therein and accompanying 
neutrino emissivities is in progress. 
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